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^ ■ Abstract 



In a contact density chain-growth study we investigate the solubihty-temperature pseudo-phase 



Q , diagram of a lattice polymer in a cavity with an attractive surface. In addition to the main phases 

^ ' of adsorbed and desorbed conformations we find numerous subphases of collapsed and expanded 
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I. INTRODUCTION 



The requirement of higher integration scales in electronic circuits, the onset of nanosen- 
sory applications in biomedicine, but also the fascinating capabilities of modern experimental 
setup with its enormous potential in polymer and surface research recently led to an increas- 
ing interest at the hybrid interface of organic and inorganic matter [l, 2, s], ^, O]. This also 
includes numerous detailed studies, e.g., of polymer film wetting phenomena 
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pat- 



, charged adsorbed 



li- 



tem —ion UU, prote,„-Mga„d b,„d,„g and docking 
polymers [13| as well as deposition and growth of polymers at surfaces 

In most theoretical and computational studies the polymer is anchored at the substrate 
with one of its ends which reduces the entropic freedom of the polymer. These surface- 
grafted polymers 1^, 16, 17, Q, 19, 20] are, e.g., of particular interest in studies of shape 
transformations j21|, e.g., as reaction to external fields 3, 23, 24] . However, in many 
recent experiments of organic-inorganic interfaces the setup is different and is more 



adequately described by a polymer moving in a cavity with one adsorbing surface [25|, |26 |. 
The main difference of such nongrafted polymers considered in this work is of entropic kind: 
In the desorbed phase the polymer can move freely within the cavity, and the polymer can 
fold into conformations, where the ends have no contact with the surface. 

This paper is organized as follows. In Sect. [Ill we describe the details of a minimalistic 
model for the hybrid system. The main result, the solubility-temperature pseudo-phase 
diagram, is presented and discussed in Sect. Illli The interpretation is consolidated by 
exemplified studies of fluctuations and correlations of relevant thermodynamic quantities 
such as numbers of contacts between monomers and monomer-substrate contacts as well as 
the gyration tensor, in the different phases. The contact numbers turn out to be adequate 
system parameters for the description of the macrostate of the system, and therefore the free 
energy in dependence of these contact numbers is subject of a detailed study in Sect. HVl 
This quantity is also useful for classifying the conformational transitions between the phases 
which are also discussed there. Eventually, we conclude in Sect. |V]with a summary of the 
main results. 
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attractive substrate 



FIG. 1: Cavity model used in this work. The lower of the two parallel surfaces is attractive to the 
polymer, the upper is steric only. The distance between the surfaces is lattice units. 

II. MINIMALISTIC MODEL FOR POLYMER ADSORPTION 



We employ a minimalistic simple-cubic (sc) lattice model 25|] which allows a systematic 
analysis of the conformational phases experienced by a nongrafted polymer in a cavity with 
one adhesive surface. An example for the cavity model is shown in Fig. [TJ The polymer can 
move between the two infinitely extended parallel planar walls, separated by a distance 
expressed in lattice units. The substrate is short-range attractive to the monomers of the 
polymer chain, while the influence of the other wall is purely steric. 

Denoting the number of nearest-neighbor, but nonadjacent monomer- monomer contacts 
by rim and the number of nearest-neighbor monomer-substrate contacts by n^, the energy 
of the hybrid system can be expressed in the simplest model as 



(1) 



where Ss and Sm are the respective contact energy scales, which are left open in the following. 
For simplicity, we perform a simple rescaling and set Ss = Eq and Em = EqS. Here we 
have introduced the overall energy scale Eq and the dimensionless reciprocal solubility s 
that controls the quality of the implicit solvent surrounding the polymer (the larger the 
s, the worse the solvent). Since contacts with the substrate usually entail a reduction of 
monomer-monomer contacts, there are two competing forces (rated against each other by 
the energy scales) affecting the formation of intrinsic and surface contacts. In this paper 
we mainly focus on the conformational transitions the polymer experiences under different 
environmental conditions. Concretely, we are interested in the dependence of energetic 



and structural quantities on temperature T and reciprocal solubility s in equilibrium. The 
probability (per unit area) for a conformation with surface and Um monomer- monomer 
contacts at temperature T and reciprocal solubihty s is given by 

where (7„^„„ = ^uaOQum + (1 ~ ^nso)9^^nm contact density and Z the partition sum. 

In this decomposition, (/^^ stands for the density of unbound conformations, whereas g^^nm 
is the density of surface and intrinsic contacts of all conformations bound to the substrate. 
Obviously, the number of the conformations without contact to the attractive substrate, 
(7"^, depends on the distance between the cavity walls. For a sufficiently large distance 
Zw from the substrate the influence of the neutral surface on the unbound polymer is small. 
For Zyj —>■ 00, however, (7"^ formally diverges. Therefore, the non-adhesive, impenetrable 
steric wall is necessary for regularization. 

We studied polymers with up to 200 monomers by applying the contact- density chain- 
growth algorithm which is an irnproved variant of the recently developed multicanonical 



chain-growth sampling method 27|, |28|. All these methods set up on a variant of the 



pruned-enriched variant 29|] of Rosenbluth sampling [301] • The main advantage of the im- 
proved method is that it directly samples the contact density gnsUm, which is very useful for 
problems, where the model provides different energy scales. This generalizes the ordinary 
multicanonical version [^^j which samples the density of states, i.e., the number of states for 
given energy. Here we can set the two independent energy scales Em and Eg or their ratio s, 
respectively, after the simulation. This allows to introduce the reciprocal solubility s as a 
second environmental parameter in addition to the temperature T. 

The partition sum of the system as a function of these two parameters is simply Z = 
J2ns,nm 9nsn„,^w{^o{^s + STim) /ksT} and the statistical average of any function Oins^rim) 
is given by the formula 

(0)(T,s)= ^ 0{ns,n^)pT,s{ns,n^), (3) 

which is very convenient since it only requires to estimate the contact density UnsUm in. the 
simulation. Denoting contact correlation matrix elements as Mxy{T,s) = {xy)c = {xy) — 
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{x){y) with x,y = ns,nm, the specific heat can be written as 



(4) 



V 



All quantities depending only on the contact numbers Um and can therefore simply be 
calculated from the estimate of the contact density provided by our simulation method. 

Although the two contact parameters are sufficient to describe the macrostate of the 
system and their fluctuations characterize the main pseudo-phase transition lines, it is often 
useful to introduce also nonenergetic quantities such as the end-to-end distance and the 
gyration tensor for gaining more detailed structural information of the polymer. For our 
specific problem at hand it is particularly useful to study the structural anisotropy of the 
adsorbed polymer in the different phases. To this end, we define the general gyration tensor 
for a polymer chain of beads with the components 

n=l 

where x\^\ i = 1,2,3, is the ith Cartesian coordinate of the nth monomer and Xi = 
J2n=i x^f^ /N is the center of mass with respect to the ith. coordinate. Anisotropy in the poly- 
mer fluctuations is connected with the system's geometry and therefore it will be sufficient to 
study the components of the gyration tensor parallel (x, y components) and perpendicular 
(in z direction) to the planar walls. 



n=l L 



(6) 



and 



^ - 2 



^1 = ^E(-^"^-^) • (7) 

n=l 

The gyration radius is then simply the trace of the gyration tensor, i?gyj, = Tri?^ = 
J2^=iRii = -Ry + R\. The calculation of statistical averages for quantities R that are not 
necessarily functions of the contact numbers Ug and rim cannot be performed via Eq. (I3l). 
In this case only the more general relation (R) = X^x -^(X) 6xp{— -E's(X)/A;5T}/Z holds, 
where the sum runs over all polymer conformations X. Introducing the accumulated den- 
sity R^c{K,n'^) = Ex^(X)5„4x)<5n„(x)n;,/^, where Sij is the Kronecker symbol, the 
expectation value can be expressed, however, in a form similar to Eq. ([3]): 

{R)= E i?accK,n^)e-^»("-"-)/'=-^. (8) 



The quantity Reccing, Um) can easily be measured in simulations with the contact density 
chain-growth algorithm. In the following we use natural units, i.e., we set /cb = = 1- 



III. THERMODYNAMIC BEHAVIOR OF THE HYBRID SYSTEM 



For our exemplified study of the hybrid system in equilibrium we chose a polymer with 
179 monomers. Since this is a prime number, the polymer is unable to form perfect cuboid 
conformations on the sc lattice, as it is, e.g., the case for a 100- mer. 25| There we found 
two low-temperature subphases dominated by the same 4x5x5 cuboid. In one subphase 
it had 20 surface contacts, while in the other the cuboid was simply rotated, entailing 25 
surface contacts. This is a typical example, where the exact number of monomers in the 
linear chain is directly connected with the occurrence of such specific pseudo-phases which 
are not, of course, phases in the traditional view. Nonetheless, the enormous progress 
in high-resolution experimental structure analyses and in the technological equipment for 
precise polymer deposition, as well as the natural finite length of classes of polymers (e.g., 
peptides and proteins), explain the growing interest in pseudo-phases and the conformational 
transitions between them. Here we mainly focus on the expected thermodynamic phase 
transitions {igI . \y\ and low-temperature higher-order layering pseudo-phase transitions jl^ . 
The following results were obtained from contact-density chain-growth simulations of the 
179- mer in a cavity with = 200 (see Fig.[l]), choosing uniformly distributed starting points 
at random. In eight independent runs 1.6 x 10^ polymer conformations were generated 
in total. The resulting contact density gnsum accumulated densities like Ra.cc{ns,nm) 
are independent of external parameters such as temperature T and reciprocal solubility s. 
Concrete values of statistical quantities for specific parameter settings are obtained by simple 
reweighting as in Eqs. ([2]) and ([S]). 



A. Solubility-temperature pseudo-phase diagram 

Discontinuities or divergences of energetic and nonenergetic fiuctuations as functions 
of external parameters reveal typically dramatic cooperative transitions in the collective, 
macroscopic behavior of the system's microscopic degrees of freedom in the thermodynamic 
limit. These transitions separate then thermodynamically stable phases and the transitions 
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FIG. 2: (Color online) Solubility-temperature pseudo-phase diagram of a 179-mer. The color codes 
the specific heat as a function of reciprocal solubility s and temperature T ~ the brighter the larger 
its value. Drawn lines emphasize the ridges of the profile and indicate transitions between the 
different conformational phases. Black lines mark expected thermodynamical phase transitions, 
while white lines belong to pseudo-transitions specific to finite-length polymers. Along the dashed 
black line coexisting desorbed and adsorbed conformations are equally probable. 



can uniquely be identified by certain values of the external parameters, e.g., the transi- 
tion temperature. Usually, all fluctuations collapse at the same parameter sets. But, this 
"traditional view" is only true in the thermodynamic limit. Finite-size systems usually 
exhibit a zoo of crossover- or pseudo-transitions, most of which disappearing in the ther- 
modynamic limit. In special cases, e.g., proteins, where the specific amino acid sequence 
is of finite length, no phase transitions in the strict sense happen at all. Still, peaks in 
curves of fluctuating quantities can be signatures for "cooperative activity" , but this is not 
necessarily indicated by all fluctuations considered, and if, then typically at different param- 
eter values. Nonetheless, in protein science, pseudo-transitions such as conformational 
transitions are important in the understanding of secondary structure formation and the 



tertiary hydrophobic-core collapse. For polymers, main 



is probably of second order, is of particular interest. 29|, |3l|] This is a real thermodynamic 



y the B collapse transition, which 



phase transition. Nonetheless, at least for finite systems, an additional first-order-like g. 



or crystallization transition at lower temperatures is also conjectured for polymers. 32 | 



assy 



All these peculiarities of finite polymer systems are also relevant for the adsorption prob- 
lem we consider here. In Fig. |2] we have plotted the projection of the specific heat profile 
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onto the solubility-temperature plane as obtained from our simulation of the 179-mer in a 
cavity with z^, = 200. The color code reflects the value of the specific heat and the brighter 
the shading, the larger the value of Cy- Black and white lines emphasize the ridges of the 
profile. Since we consider the specific heat as appropriate to identify pseudo-phases, these 
ridges mark the pseudo-phase boundaries. As expected, the pseudo-phase diagram is divided 
into two main parts, the phases of adsorption and desorption. The two desorbed pseudo- 
phases DC (desorbed-compact conformations) and DE (desorbed-expanded structures) are 
separated by the collapse transition line which corresponds to the G transition of the infinite- 



length polymer which is allowed to extend into the three spatial dimensions |33|]. The region 
of the adsorbed pseudo-phases is much more complex, and little is known about its details, 
since it is relevant at lower temperatures, where conventional Monte Carlo methods with 



pivot- 
(AE) 



ike updates usually tend to fail. The presence of general phases of adsorbed-expanded 



3J] and adsorbed-compact (ACl, AC2) conformations was postulated in adsorption 



studies of grafted polymers and the existence of an additiona 



globules (AG) 



was assumed 
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phase of surface-attached 



191]. In a recent study [191], it was argued that the 
layering transition between ACl and AC2 is a thermodynamic phase transition. Although 
the polymer in our study is still relatively small, we can clearly identify pseudo-phases in 
Fig. [2] which can be assigned these labels, too. Those regions are separated by the black 
lines indicating the transitions between them. We expect that these are transitions in the 
thermodynamic meaning; only the precise location of the transition lines will still change 
with increasing length of the polymer. Thus, this picture confirms the previously assumed 
phases and it provides evidence that the AG phase is indeed there. Furthermore, we have 
also highlighted by white lines transitions between pseudo-phases which will probably not 
survive in the thermodynamic limit. This concerns, e.g., the higher-order layering transi- 
tions among the compact pseudo-phases AC2ai^2-d. In the following sections we will analyse 
the properties of the pseudo-phases in more detail. 

B. Contact-number fluctuations 



The contact numbers and can be considered as system parameters appropriately de- 
scribing the state of the system and are therefore useful to identify the pseudo-phases. Peaks 
and dips in the external-parameter dependence of self- correlations (n^)c, {n'^)c and cross- 
correlations {nsUmjc indicate activity in the contact-number fiuctuations and, analysing the 
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FIG. 3: (Color online) Expectation values, self- and cross-correlations of the contact numbers Ug 
and Um as functions of the temperature T in comparison with the specific heat for a 179-mer in 
solvent with s = 1. 

expectation values (n^) and (rim) in these active regions of the external parameters T and s, 
allow for an interpretation of the respective conformational transitions between the pseudo- 
phases. 

In Fig. [3l we have plotted for the 179-mer these quantities and, for comparison, the 
specific heat as functions of the temperature T at a fixed solvent parameter s = 1. This 
example is quite illustrative as the system experiences several conformational transitions 
when increasing the temperature starting from T = (see Fig. [2]). At temperatures very 
close to T = (pseudo-phase ACl) all 179 monomers have contact to the substrate and 153 
monomer-monomer contacts are formed. This is the most compact contact set being possible 
for topologically two-dimensional, film-like conformations. It should be noted, however, that 
approximately 2 x 10^^ conformations (self-avoiding walks) belong to this contact set. [35| 
This high degeneracy is an artefact of the minimalistic lattice polymer model used. It is 
remarkable that the conformations with the highest number of total contacts n = ns+Um are 
film-like compact [n = 332). All other conformations we found possess less contacts, even 
the most compact contact set that dominates the five-layer pseudo-phase AC2ai (n^ = 36, 
rim = 263, i.e., n = 299). The reason is that for low temperatures, those macrostates 
are formed which are energetically favored. Entropy is not yet relevant - for the s = 1 
example (n^) drops to 149 only up to T 0.3. Increasing the temperature further, the 
situation dramatically changes, as can be seen in Fig. [3l In a highly cooperative process, 
the average number of intrinsic contacts {um) significantly increases (to ~ 208) at the 
expense of surface contacts {{ng) drops to approximately 104). Consequently, the strong 



fluctuations {nl^)c signalize a conformational transition, and the anticorrelation indicated 
by {nsnm)c confirms that surface contacts turn into intrinsic contacts, which indirectly leads 
to the conclusion that the film-like structure is given up in favor of layered, spatially three- 
dimensional conformations. The system has entered subphase AGe which is the part of 
the phase AG, where two-layer conformations dominate. The subphase transition near 
T ^ 0.7 from the two-layer (AGe) to the bulky regime of AG is due to the ongoing, rather 
unstructured expansion of the polymer into the z direction by forming so-called surface- 
attached globules 17|. This is accompanied by a further reduction of surface contacts, while 
the number of intrinsic contacts changes weakly. Approaching T 2.0, the situation is 
just vice versa. Intrinsic contacts dissolve and the system experiences a conformational 
phase transition from globular conformations in AG to random strands in AE. Crossing this 
transition line, the system enters the good-solvent regime. Eventually, close to T ^ 2.8, the 
polymer unbinds off the substrate. A clear signal is observed in the fluctuations of n^, i.e., 
the number of average surface contacts rapidly decreases. The expanded polymer is "free" 
and the influence of both walls is effectively steric. This phase (AE) is closely related to the 
typical random-coil phase of entirely free and dissolved polymers in good solvent. 

This example shows that a study of the contact number fluctuations is indeed sufficient to 
qualitatively identify and describe the conformational transitions between the pseudo-phases 
of the hybrid system. For this reason, Ug and Um are adequate system parameters playing a 
similar role as order parameters in thermodynamic phase transitions. 



C. Anisotropic behavior of gyration tensor components 

One of the most interesting structural quantities in studies of polymer phase transitions 
is the gyration tensor ([5]). For our hybrid system we expect that the respective components 
parallel ([6]) and perpendicular (I7j) to the substrate will behave differently when the polymer 
passes pseudo-transition lines. In order to prove this anisotropy explicitly, we have plotted 
in Fig. m the expectation values, and the fluctuations of these two components, 

d{R\\,i.)/dT, again for the polymer in solvent with s = 1. For interpreting the peaks of 
the fluctuations, we have also included once more the specific heat curve for comparison. 
The immediate observation is that the temperatures, where one or both gyration tensor 
components exhibit peaks, almost perfectly coincide with those of the specific heat. This is 
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FIG. 4: (Color online) Anisotropic behavior of gyration tensor components parallel and perpen- 
dicular to the substrate and their fluctuations as functions of the temperature T for a 179-mer at 
s = 1. For comparison, we have also plotted the associated specific-heat curve. 

a strong confirmation for the phase diagram in Fig. [2] which is based on the specific heat. 
Obviously, even for the rather short polymer with 179 monomers, we encounter the onset of 
fluctuation collapse near the (pseudo-) phase transitions. This is very promising for future 
quantitative flnite-size scaling analyses. 

At very low temperatures, i.e., in pseudo-phase ACl, we have argued in the previous sec- 
tion that the polymer-conformation is the most compact single-layer fllm. This is conflrmed 
by the behavior of {R\\) and the latter being zero in this phase. A simple argument 

that the structure is indeed maximally compact is as follows. It is well known that the 
most compact shape in the two-dimensional continuous space is the circle. For n monomers 
residing in it, n ~ vrr^, where r is the (dimensionless) radius of this circle. The usual squared 
gyration radius is 



R\\ ~ 5.46 of the ground-state conformation we identifled in phase ACl. Note that the most 
compact shape in the simple lattice polymer model we used in our study is a square and not 
a discretized circle. IS™ 




Near T ^ 0.3, the strong layering transition from ACl to AGe is accompanied by an 
immediate decrease of (-Ry), while {Ra_) rapidly increases from zero to about 0.5 which is 
exactly the gyration radius (perpendicular to the layers) of a two-layer system, where both 
layers cover approximately the same area. Note that the single layers are still compact. 




(9) 




Indeed, this is close to the value 




11 



but not maximally. Applying the same approximation as in Eq. ([9]), the planar gyration 
radius for each of the two layers is now (with n ^ N/2) R™'^ ^ 3.77, while we measured 
in this phase (AGe) R\\ ^ 4.05. This separates the subphase AGe from the other two- 
layer pseudo-phase AC2d in Fig. [21 where the dominating conformation has perfect two- 
layer (lattice) structure with R\\ ^ 3.85 (this is the same 2% difference between continuous 
and lattice calculation for perfect shapes as above). We will discuss the conformational 
peculiarities in the following in more detail. The subphase transition from AGe to AG near 
T ^ 0.7 is accompanied by a further decrease of (-R||) whereas {R±) increases, i.e., the 
height of the surface-adsorbed globule increases at the expense of the width. This tendency 
is stopped when approaching the transition (T ^ 2.0) from the globular regime AG to the 
phase of expanded, but still adsorbed conformations. While {R±) remains widely constant 
(the fluctuation does not signalize any transition), the polymer strongly extends in the 
directions parallel to the substrate, as indicated by the peak of d{R\\) /dT. After unbinding 
from the substrate, parallel and perpendicular gyration radii behave widely isotropically 
{{R±)'^ ~ ~ (-Rgyr)^/3) as the influence of the isotropy-disturbing walls is weak in 

this regime. 



IV. THE WHOLE PICTURE: THE FREE-ENERGY LANDSCAPE 

It was shown in Sect. IlII Bl that the contact numbers Ug and rim are unique system 
parameters for the pseudo-phase identification of the hybrid system. We define the restricted 
partition sum for a macrostate with Ug surface contacts and monomer-monomer contacts 
by 

= ^^n.n^e-^^^"-"-)/^'-^, (10) 

such that Z = Y^ns,n,n ZT,s{ns,nm)- Assuming as usual that the dominating macrostate is 
given by the minimum of the free energy as a function of appropriate system parameters, it 
is useful to define the specific contact free energy as a function of the contact numbers Ug 
and rim, 

FrA^s^nm) = -fcBTln^7„^„„e-^=("-"'")/'=^^ (11) 
= Es{ns,nm) -TS{ns,nm), 
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FIG. 5: (Color online) Map of all minima of the contact free energy Fj-^sins, rim) in the parameter 
intervals T E [0, 10] and s E [—2, 10] for the 179-mer. The solid line connects the free-energy 
minima taken by the polymer in solvent with s = 1 by increasing the temperature from T = to 
T = 5 and thus symbolizes its "path" through the free-energy landscape. The solid line is only a 
guide to the eyes. 

identifying kBliagn^nm = S{ns,nm) as a "micro-contact" entropy. For given external pa- 
rameters T and s, this relation can be used to determine the minimum of the contact free 
energy and therefore allows the identification of the dominant macrostate with respect to 
the contact numbers. In turn, this quantity allows for an alternative representation of the 
pseudo-phase diagram, complementary to the one shown in Fig. [2] in that it is related to 
the contact numbers Ug and rim- This is done by determining for (in principle) all values 
of the external parameters T and s the minima of the contact free energy (ITTl) . Then the 
pair of values Ug and ra^ of the minimum contact free-energy state are marked in an ris-rim 
phase diagram. This is shown in Fig. El where all free-energy minima of the 179-mer for 
the parameter set T G [0, 10] and s G [—2, 10] are included and, based on the arguments 
of the previous section, differently shaded according to the pseudo-phase they belong to. 
The nice thing of this representation is that it allows the differentiation of continuous and 
discontinuous pseudo-phase transitions. 

The first important observation is that the diagram is divided into two separate regions, 
the pseudo-phases of desorbed conformations (DC and DE) and the remaining different 
phases of adsorption. The "space" in between is blank, i.e., none of these (possible) confor- 
mations was found to be a free-energy minimum conformation. This shows that transitions 
between the adsorbed and desorbed pseudo-phases are always first-order-like. It should be 
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TABLE I: Representative minimum free-energy examples of conformations in the different pseudo- 
phases of a 179-mer in a cavity. The substrate is shaded in hghtgray. 



pseudo-phase 


example 




DC 


'V 


219 


DE 




50 








AG 




49 227 


AC2ai 




36 263 


AC2ao 




39 256 


AC2b 


•XBHHCI 


46 257 


AC2c 




60 251 


AC2d 




90 231 


AGe 




103 207 


ACl 




179 153 



noted, that the regime of contact pairs {ns,nm) lying above the shown compact phases is 
forbidden, i.e., conformations with such contact numbers do not exist on the sc lattice. 

The second remarkable result is that the pseudo-phases DC, DE, AE, and AG are "bulky" , 
while all AC subphases are highly localized in the plot of the free-energy minima. Comparing 
with Fig. [21 the conclusion is that conformations in the AC phases are energetically favored 
(more explicitly, for s/T > 0.8 in ACl and s/T > 2.2 in the AC2 subphases), while the 
behavior in the other pseudo-phases is entropy-dominated: The number of conformations 
with similar contact numbers in the globular or expanded regime is higher than the rather 
exceptional conformations in the compact phases, i.e., for sufficiently small s/T ratios the 
entropic effect overcompensates the energetic contribution to the free energy. 

The subphases AC2ai,2-d are strongly localized, thorn-like "peninsulas" standing out from 
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the AG regime. The discrete number and their separation leads to the conclusion that they 
have related structures. Indeed, as can be seen in Tabled where we have listed representative 
conformations for all pseudo-phases, the few conformations dominating these subphases 
exhibit compact layered structures. The most compact three-dimensional conformation with 
263 monomer-monomer contacts and 36 surface contacts is favored in subphase AC2ai and 
possesses five layers. Starting from this subphase and increasing the temperature, two things 
may happen. A rather small change is accompanied with the transition to AC2a2, where 
the number of intrinsic contacts is reduced but the global five-layer structure remains. On 
the other hand, passing the transition line towards AC2b, the monomers prefer to arrange 
in compact four-layer conformations. Advancing towards AC2d, the typical conformations 
reduce layer by layer in order to increase the number of surface contacts. In AC2d there 
are still two layers lying almost perfectly on top of each other. This is similar in subphase 
AGe, where also two-layer but less compact conformations dominate. In pseudo-phase ACl 
only the film-like surface layer remains. The reason for the differentiation of the phases 
ACl and AC2 of layered conformations is that the transition from single- to double-layer 
conformations is expected to be a real phase transition, while the transitions between the 
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higher- layer AC2 subphases are assumed to disappear in the thermodynamic limit. 

As can be seen in Fig. |2l a transition between ACl and the phase of adsorbed, expanded 
conformations, AE, is possible. Since these two phases are connected in Fig. O we expect 
that the transition in between is second-order-like. Indeed, this transition is strongly related 
with the two-dimensional G transition since, close to the transition line, all monomers form 
a planar (surface-) layer. Similarly, there is also a second-order-like transition line so(T) 
between AG and AE which separates the regions of poor (AG: s > Sq) and good (AE: 
s < So) solvent. Also, the transition between the desorbed compact (DC) and expanded (DE) 
conformations is second-order-like: This transition is strongly related with the well-known G 
transition in three dimensions. 3lf| Eventually, the transitions from the layer-phases AC2a2, 
AC2b, AC2c, and AGe to the globular pseudo-phase AG as well as transitions between 
pseudo-phases dominated by the same layer type (i.e., between the two-layer subphases 
AC2d and AGe, and between the five-layer subphases AC2ai and AC2a2) are expected to 
be continuous. 

On the other hand, the transitions among the energetically caused compact low- 
temperature pseudo-phases are rather first-order-like, due to their noticeable localization in 
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(a) 





FIG. 6: (Color online) Probability distributions PT,sins,nm) for the 179-mer in solvent with s = 1 
(a) near the layering transition from ACl to AGe at T ~ 0.34 and (b) near the adsorption- 
desorption transition from AEl to DE at T ~ 2.44. Both transitions are expected to be real phase 
transitions in the thermodynamic limit and look first-order-like. 

the map of free-energy minima (Fig. [5]). The possible transitions (see Fig. [2]) are AC2ai 2- 
AC2b, AC2b-AC2c, and AC2c-AC2d, respectively. Even more interesting, however, are the 
transitions from the single-layer pseudo-phase ACl to the double-layer subphases AC2d and 
AGe. In the previous sections we already discussed this transition for the special choice 
s = 1, where near T ^ 0.3 the fluctuations of the contact numbers and the components 
of the gyration tensor exhibit a strong activity. We have included into Fig. [5] the "path" 
of macrostates the system passes by increasing the temperature from T = to T = 5. At 
T = the system is in a film-like, single-layer state. Near T ^ 0.3 it indeed suddenly 
rearranges into two layers and enters subphase AGe in a single step. In Fig. [6t^a) we have 
plotted the probability distribution pT,s{ns, Um) for s = 1 and T = 0.34 and it can clearly be 



seen that two distinguished macrostates coexist. [37|, l38|] Increasing the temperature further, 
the system undergoes the continuous transitions from AGe via AG until it unfolds when 
entering pseudo-phase AE. The system is still in contact with the substrate. Close to a tem- 
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perature T ^ 2.4, however, the unbinding of the polymer off the substrate happens (from 
AE to DE). Comparing Figs. [5] and [6](b), where the probabihty distribution at T = 2.44 is 
shown, we see also a clear indication for a discontinuous transition. Note that we consider 
here the transition state, where the two minima of the free energy coincide 39| (see also the 
black dashed line in Fig. [2]) and not the point, where the width of the distribution, i.e., the 
specific heat, is maximal. Since the system is finite, the transition temperature (T ^ 2.8), 
as signalled by the fluctuations studied in the previous sections, deviates slightly from the 
transition-state temperature reported here. 

V. SUMMARY 

In this paper, we have studied in detail the solubility-temperature (pseudo-)phase dia- 
gram of a polymer in a cavity with an attractive substrate. We identified the thermodynamic 



phases of adsorbed compact and expanded (AC, AE) and desorbed (DC [33|, DE) confor- 
mations as well as the previously not yet clearly confirmed phase of adsorbed globules (AG). 
Although the polymer in our study possessed only = 179 monomers, these (pseudo-) phases 
are expected to be stable also in the thermodynamic limit N ^ oo. Other noticeable phase 
transitions in the compact-globular adsorbed regime (ACl-AC2d, ACl-AGe) are the ener- 
getic layering transitions from film-like surface-layer to double-layer conformations which 
are also believed to survive the thermodynamic limit. In addition, further subphases of 
higher-order layers were found in low-temperature regions and bad solvent (AC2ai 2, AC2b, 
and AC2c). The most compact three-dimensional conformation found is cube-like and forms 
five layers (in subphase AC2ai). 

The (pseudo-)phase diagram is based on the specific-heat profile as a function of temper- 
ature and reciprocal solubility. Although this profile allows for the identification of phases 
and their boundaries it does tell little about the conformational transitions between the 
phases. For this purpose we considered expectation values and fiuctuations for the numbers 
of monomer-surface contacts, n^, and intrinsic monomer-monomer contacts, Um, separately. 
These contact numbers turned out to be sufficient to describe the macrostate of the system 
and therefore are useful to describe the conformations dominating the different phases. This 
view was completed by an exemplified study of the anisotropic behavior of the gyration 
tensor components of the polymer parallel and perpendicular to the substrate. 
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Another central aspect was the classification of the conformational transitions between 
the (pseudo-)phases. Based on the contact numbers Us and rim, we defined an appropriate 
free energy and studied the distribution of the minima in the ng-rim space. Prom this kind 
of free-energy landscape, we found strong indications that the binding-unbinding transitions 
between the adsorbed and desorbed phases are first-order-like. This was also observed for 
the layering transitions. On the other hand, the transitions across the line separating good 
and poor solvent, i.e., between the compact (or globular) and the expanded conformations, 
are rather second-order-like. This is in coincidence with the known behavior of free polymers 
at the © collapse transition in two and three dimensions. 

Since the experimental equipment and the technological capabilities have nowadays 
reached an enormous standard of high single- molecular resolution, we expect that it should 
be possible to verify experimentally not only the existence of the described thermodynamic 
phases, but also the pseudo-phases being only relevant for finite polymers and specific to 
their lengths. 
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